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Abstract 

Some extensions of the Standard Model provide Dark Matter candidate par- 
ticles which can have a dominant coupling with the lepton sector of the ordinary 
matter. Thus, such Dark Matter candidate particles (x°) can be directly de- 
tected only through their interaction with electrons in the detectors of a suitable 
experiment, while they are lost by experiments based on the rejection of the elec- 
tromagnetic component of the experimental counting rate. These candidates can 
also offer a possible source of the 511 keV photons observed from the galactic 
bulge. In this paper this scenario is investigated. Some theoretical arguments are 
developed and related phenomenological aspects are discussed. Allowed intervals 
and regions for the characteristic phenomenological parameters of the considered 
model and of the possible mediator of the interaction are also derived considering 
the DAMA/Nal data. 
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1 Introduction 

Dark Matter particles with dominant interaction on electrons have been considered 
in literature [1, 2, 3, 4]. In particular, from a phenomenological point of view, Dark 
Matter (DM) candidates with electron interactions can offer possible sources for the 
511 keV positron annihilation line observed from the galactic bulge [5, 6]. These 
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candidates can be either light (MeV scale) [1] or heavy (GeV or larger scale) [2, 3]. 
They are expected to interact with electrons both through neutral light (MeV scale) U 
or Z' bosons or through heavy charged mediators (which can eventually be nearly 
degenerate with x°) [3]. Recently data collected by some accelerator experiments 
have been analyzed in terms of a ~ 200 MeV neutral boson which couples to quarks 
with flavour changing transition: s — > d/i + ji~ [7, 8]. Other results showing some 
resonances at energies lower than the two-muon [7] and the two-pion [9] disintegration 
thresholds have been associated with a Goldstone neutral boson of ~ 20 MeV mass. 
Moreover, some excess has been achieved in dedicated experiments on low energy 
nuclear reactions searching for possible e + — er pairs driven by the presence of a 
neutral boson with a mass around 10 MeV [10]. 

Let us remark that in the frameworks where the mediator is cither ail charged 
boson or a neutral boson providing a flavour changing transition among quarks - the 
elastic scatterings of the DM candidate \° particles on nuclei would be either forbidden 
or suppressed; hence, the scattering on electrons would remain the unique possibility 
for the direct detection of the x° particles. 

On the other hand, from a pure theoretical point of view, it is also conceivable 
that the mediator of the DM particle interactions can be coupled only to the lepton 
sector of the ordinary matter. Thus, in this case the DM particles can just interact 
with electrons and cannot with nuclei. This is suggested in ref. [4] for the U boson 
and can also be the case of some extensions 2 of the Standard Model providing a 
quark-lepton discrete symmetry SU(3)i x SU(3) q x SU(2)l x U(l). In these latter 
models, leptons (as well as quarks) are assumed to have three "leptonic (I) colours" 
and to interact through the gauge group SU (3);, analogously as the QCD colour group 
SU(3) g . Moreover, at some high energy scale a symmetry breaking SU(3)i — > SU(2)' 
is expected, giving high mass to the "exotic" leptonic degree of freedom and leaving 
light the "standard" leptons [13]. In these scenarios, the heavy "exotic" leptonic 
degree of freedom provides both heavy charged ±1/2 fcrmions, which are expected to 
be confined into exotic leptonic hadrons by the unbroken gauge group SU{2)' [13], and 
heavy neutrinos [12, 13]; hence, they can be considered as Dark Matter candidates 
with dominant interaction on electrons. 

Moreover, it is worth to note that other possibilities can exist. For example, super- 
symmetric (SUSY) theories can offer configurations in the general SUSY parameter 
space where the lightest supersymmetric particle (LSP) has an interaction with elec- 
tron dominant with respect to that with quark. 

These DM candidate particles can be directly detected only through their interac- 
tion with electrons in the detectors of a suitable experiment, while they are lost by 
experiments based on the rejection of the electromagnetic component of the experi- 
mental counting rate. 

In the present paper this kind of DM candidates are investigated, some theoretical 
arguments are developed and related phenomenological aspects are discussed. In par- 
ticular, the impact of these DM candidates will also be discussed in a phenomenological 
framework on the basis of the 6.3 a C.L. DAMA/Nal model independent evidence for 
particle Dark Matter in the galactic halo [14, 15]. We remind that various corollary 

2 For example from the extended Pati-Salam gauge group SU(6) X SU{2)l x SU(2)r [11] or from 
[St/(3)] 4 quartification [12]. 
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analyses, considering some of the many possible astrophysical, nuclear and particle 
Physics scenarios, have been analysed by DAMA itself both for some WIMP/WIMP- 
like candidates and for light bosons [14, 15, 16, 17, 18, 19], while several others are 
also available in literature, such as e.g. rcfs. [20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. 
Many other scenarios can be considered as well. At present the new second generation 
DAM A/LIBRA set-up is running at the Gran Sasso Laboratory. 

2 Detectable energy in x° ~ electron elastic scatter- 
ing 

The practical possibility to detect electron interacting DM candidates (hereafter x° 
with mass m x o and 4-momentum k^) is based on the detectability of the energy released 
in x° - electron elastic scattering processes (see Fig. 1). 




Figure 1: The x° ~ e~ elastic scattering and definition of the momentum variables in 
the laboratory frame. In the text a contact interaction has been assumed (also see 
Appendix B) as suitable approximation of the process. 

Generally, these processes are not taken into account in the DM field since the 
electron is assumed at rest and, therefore, considering the \° particle velocity \v x o\ ~ 
300 km/ s, the released energy is of the order of few eV, well below the detectable energy 
in any considered detector in the field. However, the electron is bound in the atom 
and, even if the atom is at rest, the electron can have not negligible momentum, p. For 
example, the bound electrons in Nal(Tl) offer a probability equal to <~ 1.5 x 10~ 4 to 
havep > 0.5 MeV/c; such a probability is quite small, but not zero. Hence, interactions 
of x° particles with these high-momentum electrons in an atom at rest can give rise 
to detectable signals in suitable detectors. In particular, after the interaction the final 
state can have - beyond the scattered x° particle either a prompt electron and an 
ionized atom or an excited atom plus possible X-rays/ Auger electrons. Therefore, 
the process produces X-rays and electrons of relatively low energy, which arc mostly 
contained with efficiency ~ 1 in a detector of a suitable size. Thus, the total detected 
energy, Ed = fco — k' = p' a — po (where fco, k' , p' and pg are the time components 
of the respective 4- vectors in the laboratory frame, see Fig. 1), can be evaluated 
considering the energy conservation in the centre of mass (CM) frame of the \° ~ e ~ 
system. Defining j3 = , fc ^ p as the velocity of the CM frame with the respect to the 
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laboratory frame and 7=1/ \J\ — (3 2 Lorentz boost factor, one can write the energies 
of the electron before and after the scattering by using the variables in the CM frame 
through the Lorentz transformations: 



Po = l(Po,CM + P ■ Pcm) and p' = 7(Po,CM + ■ p' CM ). 



(1) 



Since we are dealing with elastic scattering, po,cm = Pqcm an d |pcm| = Ip'cmL so 
that, by subtraction, one obtains: 



E d = 7 ( /3 • P 'cm ~ ■ Pcm) =7/3 Pcm (cos9' - cosO) 



(2) 



where 9' is the angle between (3 and p 1 cm-, ® ^ s the angle between (3 and pcm and 
Pcm = j(jp- 0Po)- 

Therefore, fixing the input momenta of the \° particle (k) and of the electron (p), 
the maximum detected energy is given by: E+ = 7 (3 pcm (1 — cos9). Few examples of 
the dependence of E+ on the x° mass are given in Fig. 2 as function of the electron's 
momentum and of the \° velocities for head-on collisions (9 = ir). The Fig. 2 also 
points out that \° particles with m x a larger than few GeV can provide sufficient energy 
to be detected in a suitable detector. 



> 



Li 



p = 5 MeV/c 








p = 1 MeV/c 



p = 0.1 MeV/c 



m xo (GeV) 

Figure 2: Few examples of the dependence of the maximum released energy, on 
the x° mass for electron's momenta of 0.1, 1 and 5 MeV/c, for v x o ranging in the 
interval 1 2 x 10~ 3 c and for head-on collisions (6 — ir). 



It is interesting to explore two limit cases (remind that owing to the typical x° 
velocities, ko ~ m x o and k ~ m % o ■ v x o; hereafter c = 1): 

a) p <C (3m e ~ keV, that is target nearly at rest 3 : E + ~ 2j3 2 m e ~ eV. 

3 We note that in general for a target of mass my nearly at rest: E+ ~ 2/3 2 my = 

1 9 4m q iwy 

jm„ot) „ ■ -, —, r*. that is one gets the formula describing for example the WIMP-nucleus elastic 

2 x x" ( m x o+ m T)- s 

scattering. 
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b) k~3> p /3po <~ kcV; in this case one obtains pcm — P , P — v x o and, therefore, 
9 is also the angle between p and k. Hence: E + ~ v x op(l — cosO). This is the 
case of interest for the direct detection; in fact, for m x o larger than few GeV k 
is larger than the maximum momentum of a bound electron in the atom due to 
the finite size of the nucleus (~ 15 MeV in Iodine). 

In conclusion, x° particles with mass > few GeV, interacting on bound electrons 
with momentum up to ~ few MeV/c (see case b), can provide signals in the keV energy 
region detectable by low background and low energy threshold detectors, such as those 
of DAMA/Nal (see later). 



3 Cross section and counting rate 

3.1 The cross section at fixed electron momentum 

The differential cross section for x° - electron elastic scattering can be written as: 



\ M \ 1 ,„ n4c4,, ,/ ^ d 3 p' d 3 k' 

do = - — (2tt 4 (P (k + p - k' - p V r£- , j- r^— . (3 



There |M| 2 is the averaged squared matrix element and f( x o e ) is the relative velocity 
between x° and the electron. 

Integrating over d 3 k' and over the p' solid angle and considering that p' dp' = p' dp' = 
p' Q dEd, one can write: 



±L = ^Ml J-. 9{ E + -E A ). (4) 

dE d 327rw( x o e )A:oPo \k + p\ 

The Heaviside theta function defines the domain of the differential cross section. 

It is useful in the following to define the x° cross section on the electron at rest 
(p = 0); thus, one can write: 



da 
dE~ d 



|A/h( " " ! 1 0(E + - E d ) = ^6{E + - E d ), (5) 



(p 



=0 j 32nv x ok m e k E + 



| M |2 

where E + (p = 0) = 2m e v^. ~ eV and a® — 167T ^ l 2° > ■ In the following, for simplicity, 
we define a e = , then a e (p = 0) = cr°. 

3.2 The cross section for atomic electrons 

Let us now introduce in the previous evaluations the momentum distribution of the 
electrons in the atom, p(p) (see Appendix A). In particular, from eq. (4) - that is for 
a fixed p value - one can write for the atomic case: 



dE d 32nv( x o e) k p \k + p\ 
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Introducing the a e definition and replacing E+ with its expression, it is possible to 
write for the relevant case of direct detection (k 3> p 3> m e v x a): 



4^r - ° eV p(p)d<j)dcos6 e[v x op(l - cos9) - E d ]dp; (7) 

here the polar axis has been chosen in the direction of k. 



The integration over (p simply gives 2ir considering that |M| 2 does not depend on 
and that atoms with full shells (as Na + and I~) have isotropic distributions p{p)- 



3.3 The counting rate 

The expected interaction rate of x° particle impinging on the electrons of an atom can 
be derived as: 

dR p x o f da . ,o 

dE d = ^J dE d V ^ V ^ dv ^ < 8 > 

where: i) p x o ~ £,Po with p local halo density and £ < 1 fractional amount of x° 
density in the halo; ii) f(v x o) is the x° velocity (v x o) distribution in the Earth frame; 
iii) rj e is the electron's number density in the target material. 

In the reasonable hypothesis that a e does not depend on cos9, the integrand in eq. 
(8) can be evaluated considering that: 

da 2ira e p 2 ... . . . . . 

— - • v {x a e) = —g P(P){v x o - v min )0(v x a - v min )dp , (9) 

ah,d v x<) po 

where v m i n — is the minimal x° particle velocity in order to provide an energy E d 
released in the detector. 

The matrix element \M\ 2 - as well as a e in eq. (9) - can generally depend on p 
and v x o. Thus, in order to evaluate it, it is necessary to consider a specific particle 
interaction model (see Appendix B). 

For simplicity we will consider a 4-fermion contact interaction (e.g. a mediator 
with mass larger than many MeV, neglecting the 4-momentum transferred into the 
propagator). Thus, for the cases of pure V ± A and pure scalar interactions - which 

2 

are addressed in the following - one gets: a e ~ a e^- Other interaction models 
are possible and can be investigated in the future. It is worthwhile to stress that - 
although the calculations are made for the V ± A and for the scalar 4-fermion contact 
interactions - same results can be achieved for any kind of DM candidate interacting 
with electrons and with cross section a e having a weak dependence on p and v x o , that 
is a e <~ (7°. 

Finally, the expected interaction rate can be written as: 
dR £a° 2np 



dE d m x o m 



2 



Ve P POP(P) ■ I(Vmin) dp, (10) 



G 



where - pointing out the time dependence of f(v x o) - we have introduced the useful 
function: 



^—(v x o-V min )d 3 V x o ~ Io(Vmin)+Im(Vmin)-COSw(t—to).(ll) 



Here roughly to ~ 2 nd June and u> — with T = 1 yr. The cut-ofp of the halo 
escaping velocity is included into the f(v x o) function distribution. 

Therefore, the expected counting rate accounting for the energy resolution of the 
detector can be written as: 



COSU>(t — to) 



(12) 



where So and S m are the unmodulated and the modulated part of the expected signal, 
respectively. The G(E,Ed) kernel generally has a gaussian behaviour. 

Finally, we note that - since m x a is larger than few GeV (so that k ^> p) - the 
expected counting rate has a simple dependence upon a® and m x o ; therefore, the ratio 

-^—^ is a normalization factor of the expected energy distribution. 

The momentum distribution of the electrons in Nal(Tl), p(p), has been depicted 
in Fig. 3a); it has been calculated from the corresponding Compton profile, J(p), 
reported in ref. [32]. For this purpose, due to the isotropic distributions of Na + and 
I (ions with full shells) the relation J(p) = 2n f p °° p{q)qdq has been used [33, 34]. At 
high momentum the p{p) function follows the hydrogenic behaviour of the Is internal 
shell of the Iodine atom: p(p) oc (p] + p 2 ) with pj ~ 200 keV. 
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Figure 3: a) Behaviours of p(p) (solid black line) for Nal(Tl) and Io and I m for = 3 
keV in the considered halo model, A5 of ref. [31, 14]; see also text. The functions 
Io and I m are in arbitrary units, b) Behaviours of p 2 pop(p)I m for Nal(Tl) at three 
different values of the released energy: Ed — 3, 6 and 12 keV in the considered halo 
model, A5 of ref. [31, 14]; they show as the main contribution to the counting rate in 
Nal(Tl) detectors with energy threshold at 2 keV comes from electrons with momenta 
around few MeV/c. 
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As an example, in Fig. 3a) the behaviours of Io{vmin), Im(v m in) and p(p) are 
compared as function of the electron's momentum, p, for Nal(Tl) as target material 
and for the given released energy: Ed — 3 keV. In this figure as template the considered 
halo model is the A5 model of ref. [31, 14], that is a NFW halo model with local velocity 
equal to 220 km/s and density equal to the maximum value (po = 0.74 GeV cm -3 ). 




2 4 6 8 10 12 14 16 18 20 

E (keV) 



Figure 4: An example of the shapes of expected energy distributions in Nal(Tl) due 
to \° interactions with electrons for the scenario given in the text; the solid line gives 
the behaviour of the unmodulated part of the expected signal, So, while the dashed 
line is the behaviour of the modulated part, S m . In this example the normalization 
factor is J^- = 7x 10~ 3 pb/GeV. The vertical line indicates the energy threshold of 
the DAMA/Nal experiment. 

It is possible to see that - due to the behaviour of the momentum distribution of the 
electrons, p(p), at highp and due to the behaviour of the / function at low p (related to 
the f(v x o) behaviour at high velocity) - the main contribution to the counting rate in 
Nal(Tl) detectors with energy threshold at 2 keV comes from electrons with momenta 
around few MeV/c (see Fig. 3b). It is worthwhile to note that similar behaviours can 
also be obtained by using other choices of the halo model. 

Finally, an example of the shapes of expected energy distributions in Nal(Tl) due 
to x° interactions with electrons for the A5 halo model (a NFW halo model with local 
velocity equal to 220 km/s and density equal to the maximum value, see ref. [31, 14]) 

is reported in Fig. 4. In this example the normalization factor is = 7 x 10 

pb/GeV. 

4 Data analysis and results for electron interacting 
DM candidate in DAMA/Nal 

The 6.3 a C.L. model independent evidence for Dark Matter particles in the galactic 
halo achieved over seven annual cycles by DAMA/Nal [14, 15] (total exposure ~ 
1.1 x 10 5 kg x days) can also be investigated for the case of an electron interacting DM 
candidate (in addition to the other corollary quests already mentioned in the previous 
footnote 4). 
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In the analysis presented here, the same dark halo models and related parameters 
given in table VI of ref. [14] have been used; the related DM density is given in table 
VII of the same reference. Moreover, here r\ e — 2.6 x 10 26 kg -1 and the halo escaping 
velocity has been taken equal to 650 km/s. 

The results are calculated by taking into account the time and energy behaviours of 
the single-hit experimental data through the standard maximum likelihood method 4 . 
In particular, they are presented in terms of the allowed interval of the ^jj- 6 - parameter, 

obtained as superposition of the configurations corresponding to likelihood function 
values distant more than 4a from the null hypothesis (absence of modulation) in each 
one of the several (but still a very limited number) of the considered model frameworks. 
This allows us to account for at least some of the existing theoretical and experimental 
uncertainties (see e.g. in ref. [14, 15, 16, 17, 18, 19] and in literature). 

For these scenarios the DAMA/Nal annual modulation data gives for the considered 
X° candidate: 1.1 x 1CT 3 pb/GeV < < 42.7 x 1(T 3 pb/GeV at 4ct from null 

hypothesis. In particular, Fig. 5 shows the DAMA/Nal region allowed in the (£<7g vs 
m x o) plane for the same dark halo models and related parameters described in ref. 
[14]. ' 




500 



1000 



1500 



2000 



m xo (GeV) 



Figure 5: The DAMA/Nal region allowed in the (£<7g vs m x o) plane for the same 
dark halo models and related parameters described in ref. [14]. The region encloses 
configurations corresponding to likelihood function values distant more than 4cr from 
the null hypothesis (absence of modulation). We note that, although the mass region 
in the plot is up to 2 TeV, x° particles with larger masses are also allowed. 



We would like to stress that - although the above mentioned calculations have been 



4 Shortty, the likelihood function is: L = Ily^e **<i* j$* fc | > wnere ^ijh i s the number of events 
collected in the i-th time interval, by the j-th detector and in the fc-th energy bin. 2V«jfc follows a 
Poissonian distribution with expectation value Hijk = [bjk + So.fc + >Sm,k ' cosu)(ti — to)]Mj AtjABe^j,. 
The unmodulated and modulated parts of the signal, Sq % an d S m kcosu)(ti—to), respectively, are here 

functions of the only free parameter of the fit: the ratio. The bjk is the background contribution; 

Atj is the detector running time during the i-th time interval; ejj, is the overall efficiency and Mj is 
the detector mass. 
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made for the V±A and for the scalar 4-fermion contact interactions - the results given 
here hold for every kind of DM candidate interacting with electrons and with cross 
section a e having a weak dependence on p and v x o, that is a e ~ cr"; in such a case, the 

DAMA/Nal annual modulation data gives: 1.6 x 1(T 3 pb/GeV < < 53.4 x 10~ 3 

pb/GeV at 4cr from null hypothesis. 

Let us now comment some phenomenological implications about the possible me- 
diator of the interaction (hereafter U boson). The hypothesis of 4-fermion contact 
interaction still holds for U boson masses, Mjj, larger than the transferred momentum 
(Mjj > 10 MeV). In the pure V±A and pure scalar scenario, the cross section is given 
by (see Appendix B): 

o _ gfg l6G 2 m 2 x0 m 2 e _ G 2 m 2 e _ gc^m* 
° e ~ 16nm 2 x o ~ l6nm 2 x o tt ~ -kM^ ' ^ ' 

The effective coupling constant, G, depends on the couplings, c e and c x o , of the U boson 
with the electron and the x° particle, respectively. We note that limits on c e have been 
achieved by the experimental constraints on the possible U boson coupling to electron 
arising from the g e — 2 measurements: c e < 10~ 4 [4]. Moreover, more restrictive 
limits have been obtained under the assumption of universality (c M ~ c e ~ c v ) by 
considering the — 2 and v — e scattering data: < 3 x 10~ 6 ^ [4]. 

The DAMA/Nal allowed region of Fig. 5 requires values of c e well in agreement 
with these experimental upper limits. In fact, from Fig. 5 and reminding that £ < 1 
and m x o > few GeV (see above), we obtain that ct° > 10~ 2 pb. Requiring that the 
theory remains perturbative (that is, c x o < ^/in) and for Mjj ~ 10 MeV, the values 
of c e allowed by DAMA/Nal data are (see eq. (13)): c e > 5 x 10~ 7 , in agreement with 
the experimental upper limits. 




500 1000 1500 2000 



m xo (GeV) 

Figure 6: Region of U boson mass allowed by present analysis and by the g e — 2 
constrain [4] considering that £ < 1 and that the theory is perturbative (c x o < y/iir). 
See text. There U boson with Mjj masses in the sub- GeV range required by the 
analyses of ref. [1, 4, 7, 8] is well allowed for a large interval of m x o. 
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More in general, considering the limit on c e from g e — 2 data and the obtained lower 
bound ^ > 1.1 x 1(T 3 pb/GcV from the DAMA/Nal data, the allowed U boson 

masses are: Mu(GeV) < m p™ -^, as reported in Fig. 6. There U boson with Mjj 

masses in the sub-GeV range required by the analyses of ref. [1, 4, 7, 8, 9, 10] is well 
allowed for a large interval of m x o . 

5 Conclusions 

In this paper, the scenario of a DM particle x° with dominant interaction with elec- 
trons has been investigated. This candidate can be directly detected only through 
its interaction with electrons in suitable detectors. Theoretical arguments have been 
developed and related phcnomcnological aspects have been discussed. In particular, 
the impact of these DM candidates has also been analysed in a phenomenological 
framework on the basis of the DAMA/Nal data. 

For the considered dark halo models the DAMA/Nal data support for the x° can- 
didate: 1.1 x 10" 3 pb/GcV < J^- < 42.7 x 1(T 3 pb/GeV at 4cr from null hypothesis. 
Allowed regions for the characteristic phenomenological parameters of the model have 
been presented. The obtained allowed interval for the mass of the possible mediator of 
the interaction is well in agreement with the typical requirements of the phenomeno- 
logical analyses available in literature. 

Finally, we further remind that the U boson interpretation is not the unique one 
since, for example, there are domains in general SUSY parameter space where LSP- 
electron interaction can dominate LSP-quark one. 



APPENDIX 



A x interaction with atoms 

The inclusive scattering of x° particle on an atom A is here analyzed: x° A — > x°A, 
where X denotes the final state of the atom. The cross section of the process is 
obtained by summing over the possible contributions of all the X final states: 

da x o A cx^|r^| 2 ^^(Ax (A:)|x (fc')^>(^X (fc')lx W,yl) ; (14) 
x x 

here Tax is the transition amplitude when the final state is X. 

Since it has been assumed that the interaction of x° with the electrons is dominant, 
we can use a full set of electronic plane wavefunctions, e(p), and rewrite: 

(A, X °(k)\=J2(A\e(p))(e(p),x°(k)\ (15) 
p 

\x°(k'),X) = ^(e(p')\X)\ X °(k'),e(p')) . (16) 
p' 
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Therefore: 

T AX = J2(A\e(p))T(P+k- P '-k'M P ')\X) (17) 
p,p' 

where T( p+fc _p/_ fe /) = (e(p),x (k)\x°(k'),e(p')) oc M x5(p+k—p' — k') is the transition 
amplitude for free electron — x° elastic scattering and M is the matrix element reported 
in eq. (3). 

Since X is whatever final state: ^2 X (e(p')\X)(X\e(p")) = 5(p' — p"); therefore, eq. 
(14) can be written as: 

J2 T ax = ^ e {p)) T (p+k-p'-k')Tl p n, +k _ p ,_ k , ) {e{p'')\A) 

x P,p',p"' 

oc ^2p(p)\M\ 2 S(p + k-p' -k') (18) 
p,p' 

where p(p) = \ (A\e(p))\ 2 is the momentum distribution function of the electrons in the 
atom A. Finally, we can deduce da x o A — da x o e ■ p(p)d 3 p, where da x o e is the x° ~ e ~ 
elastic scattering cross section given in eq. (3). 



B The invariant amplitude for x° ~ e elastic scat- 
tering 

In the following we consider the elastic scattering of the x° fcrmion on electron by 
using a Fermi-like 4-fermion contact interaction. 

B.l The VA subcase 

The squared matrix clement, averaged over the initial spins and summed over the final 
ones, can be written as: 

\M^\* = Cl 2 L% 0) L$, (19) 

where: 

L (x°) = 111 [U x o(k')r(g v + g Al 5 )U x o(k)] [U x o(k)Y(gv + g A J 5 )U x o(k')](20) 

spin 



L tl = \ E [Ue(p'h„(c V + CAl 5 )U e (p)] [Ue{ P )l„{cv + C A1 5 )U e (p')} . (21) 
spin 

Let us focus just on eq. (20), since eq. (21) has the same structure. One can write: 



L^o, = lTr[(j(' + m x o)r(gv + gA7 5 )()i + m x o)^(g v +g Al 5 )] 



rpAA _|_ rpVA _|_ rpAV _|_ rpVV j"22) 
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-vv 



The four terms can be explicited as: 

T AA = \Tr [( ft' + m x o ) 7 <W ^ + m x o ) 7 W] (23) 

l -Tr [(#' + m x o)^g v m + m x o)^g v ] (24) 

T AV - l -Tr [(#' + m x0 )^g A fm + m x o)^g v ] (25) 

T VA = l -Tr [(#' + m x o)rg v (ji + m^g^ 5 ] (26) 
By using trace theorems one gets: 

T AA = l -g A Tr [tf-f ^ - mj^'V] = 2g\{k'»k v +k w k»-k'kg» v -m\ a g» v ){27) 

T vv = \glTr [fH'f $Y +m^Y] = 2g v (k^k v +k' v k^-k'kg^+m 2 x og^){28) 

T AV = l -Tr [(#' + m x o)r((i ~ m x o)Y W '9v] (29) 

T VA = \gv9ATr [ 7 5 li'-f /fry"] = T AV (30) 

T VA + T AV = -gyg^i^k^kp (31) 

Thus, one can write: 

L^ 0) = 2{g 2 v +g 2 A )[k'^k v + k' v k^-k l kg^} + 

+2{g 2 v - g 2 A )m 2 x0 g^ - Ag v g A ie a ^k' a k p (32) 
Finally, the matrix element for the process can be written as: 



\M VA \> = 8G 2 [A(p'k')(pk) + B{p'k){pk') - C{kk')ml - D( PP ')m 2 x o] , (33) 
where: 

A = (9v + 9a)( c v + c a) + ^9v9aCvCa = (cv9v + c A g A ) 2 + {c v g A + c A g v ) 2 
B = (.9v + 9a)( c v + c a) - ^9v9acvc a = (c v g v - c A g A ) 2 + (c v g A - c A g v ) 2 
C = (9 2 v + 9 2 a)(c 2 v ~c 2 a ) 

D = (9v - 9a)(4 + ca) (34) 
In the case of V ± A interaction (|cy| = \c A \ and \gy\ = \g A \) the matrix element 

is: 

\M v±A \i = 8G 2 [A(p'k')(pk) + B(p'k) (pk 1 )] (35) 

knowing that %° is not relativistic (see text), one obtains: (p'k')(pk) ~ p' k' poko and 
(p'k)(pk') ~ p' k' p k ; moreover for £^ <~ fceV one has p ~ po ; giving: 



\M v±A \i ~ 16G£ ±A ro£opg, ( 36 ) 



x" 

where the Fermi effective coupling constant is: Gy ±A = G 2 (cy + c A )(g v + g\). For 
this particular case, the dependence on v x a can be neglected, while the dependence on 
p are included in: p^ = p 2 + m 2 . 
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B.2 The SP subcase 

Similarly as above, one has: 

|Ms7F = G 2 L (x0) L (e) (37) 
L (x°) = \ E [U x o(k')(g s + zg Pl 5 )U xa (k)] [U x o(k)(gs+ig P j 5 )U x o(k')] (38) 

spin 



L( x o) = ^Tr [(ft' + m x a)(g s + ig P j 5 )(k' + m x a)(g s + igpj 5 )] 

nSS , rpSP , rpPS , rpPP 



T i, b +Tb P +T P b +T PP ^ 



There: 



T SS = \g 2 s Tr [(ft' + m x o)(ft + m x o)] = 2g%(k'k + m 2 x „) (40) 

T PP = -\gpTr [(ft' + m x «h 5 (ft + m x o) 7 5 ] = 2g 2 p (k'k m^) (41) 

T PS = \ig P g S Tr [(ft' + m x0 )f(ft + m x o)] (42) 

t sp + t ps = igpgsTr + mx o) 7 5 m x o] = (43) 
Hence: 

L {x0) =2[(g%+gl)k'k + (g%-g%)m x0 ] = 2(g+k'k + g-m 2 x o), (44) 
where g+ = g 2 s + g% > and = g% - g%. Finally: 



\M S p\ 2 = 4G 2 [g+c + (k'k)(p'p) + g+c-(k'k)m 2 e + g-c + (p'p)m 2 x0 + g-C-m^m 2 ] (45) 
In the particular pure scalar case (gp = cp = 0) one obtains: 



\M S \ 2 = AG 2 g 2 c 2 [(k'k) + m 2 x0 ] [(p'p) + m 2 ] ~ 8G|m % o [p' Po - p'p + m 2 e \ (46) 

Thus, considering the momentum distribution of atomic electron, for Ed <~ keV prac- 
tically p' <~ — p and, therefore: 



|M S | 2 ~ 16G 2 s m x0 p 2 (47) 

where the Fermi effective coupling constant is: G 2 S = G 2 <? s g 2 s . 

Also in this case there is a negligible dependence from v x a and a weak dependence 
from p. 
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